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Abstract 



We study scaling properties and topological aspects of the 2-d 0(3) non-linear cr-model 
on the lattice with the parametrized fixed point action recently proposed by P. Hasenfratz 
and F. Niedermayer. The behavior of the mass gap confirms the good properties of scaling 
of the fixed point action. Concerning the topology, lattice classical solutions are proved to be 
very stable under local minimization of the action; this outcome ensures the reliability of the 
cooling method for the computation of the topological susceptibility, which indeed reproduces 
the results of the field theoretical approach. Disagreement is instead observed with a different 
approach in which the fixed point topological charge operator is used: we argue that the 
discrepancy is related to the ultraviolet dominated nature of the model. 



'^Work supported in part by Fondazione "A. Delia Riccia" (Italy). 

2 On leave from the Department of Physics, University of Pisa and I.N.F.N., Pisa. 



1 Introduction 



Simulating a theory on the lattice is the most (and often the only) reliable mean to get an insight 
into its non-perturbative aspects. Once the theory has been put on the lattice, the non-trivial 
question is how to correctly reach the continuum limit. The standard procedure of adjusting the 
bare parameters of the lattice theory closer and closer to the critical values (where the lattice 
theory is expected to reproduce the continuum) is plagued by critical slowing-down and finite-size 
effects: in practice one is forced to simulate the theory at finite cut-off, where the systematic 
effects of discretization, heavily coming into play, can be removed only by an - often ambiguous - 
extrapolation. A way to reduce lattice artifacts is to use perturbatively improved lattice actions Q] , 
whose performance in numerical simulation is not, however, under theoretical control. 

A radical solution to the problem j!| is to use perfect actions, i.e. lattice actions whose spec- 
trum is completely free of lattice artifacts. Wilson's renormalization group theory || ensures that 
the fixed point (FP) of a renormalization group (RG) transformation and all lattice actions of the 
renormalized trajectory (RT) - which is the asymptotic flux line under repeated RG transforma- 
tions - are perfect (quantum) actions in the above mentioned sense. The FP action, in particular, 
represents the classical perfect action, having the same classical properties of the continuum the- 
ory §. 

A method for the determination of the FP of asymptotically free theories has been proposed 
in a recent paper j2j. The procedure has been applied to the 2-d 0(3) non-linear cr-model on the 
lattice, and a parametrization of the FP action App suitable for numerical simulations has been 
found. Strong numerical evidences [|| indicate that the action App is a good approximation of the 
perfect action at small correlation lengths also. This agrees with the general remark by Wilson Q 
that the FP action is (quantum) perfect at f-loop order in perturbation theory. A formal argument 
for this statement is presented in (j| and a further support comes from a specific 1-loop calculation 
in the 0(3) cr-model [§. 

In the first part of this paper, we test the scaling properties of the FP action of Ref. up to 
very small correlation lengths: we study the physical scaling of the mass gap - which is the only 
relevant spectral property in the 0(3) <7-model - by comparing different lattice definitions; we also 
check the rotation invariance. 

The second part of this work is devoted to the study of the topological properties of the model 
with the FP action: in particular, we address ourselves to the problem of extraction from the lattice 
of the topological susceptibility of the O (3) cr-model @@§@0©@@- We present different 
approaches to the subject - the field theoretical jl5|], the geometrical J?j and the cooling [p||l^] 
methods - and compare their outcomes. In this respect, we argue that the cooling method can be 
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safely applied with the FP action, since it possesses scale invariant classical solutions and no loss 
of topological signal is expected during the procedure of minimization. The scaling of topological 
susceptibility is a rather involved matter since the semiclassical approximation [p~T| |l8|| and some 
numerical evidences |k| [^0| indicate that the topology of the model is UV dominated. 



2 The model 



The 2-d O(N) cr-models are 0(N)-symmetric renormalizable quantum spin field theories, asymp- 
totically free for N> 3. 

They are described by the Lagrangian: 

where 4>(x) is a N-component real field satisfying the constraint <j> ■ <j> = 1 and /3 is the inverse of 
the coupling constant g. 

The infrared charge singularity of these models is responsible for the disintegration of the 
Goldstone vacuum ]2l|| ; for N> 3 the true vacuum is O(N) symmetric and non-degenerate. In the 
limit N— > oo the model contains an isovector N-plet of free massive particles only [^2)|2^]. The 
interaction of these particles is of order 1/N and so for N sufficiently large no bound states are 
present. Strong theoretical evidences bear the conjecture that the situation is the same at all N 
> 3 j24|. Consequently, the spectrum of the 0(3) cr-model is expected to consist of a single triplet 
of particles. 

Among the 0(N) cr-models, the 0(3) model in particular plays an important role in quantum 
field theory because it resembles, besides asymptotic freedom and spontaneous mass generation, 
another aspect of the 4-d non-Abelian gauge theories, i.e. the non-trivial topology. 

The topological charge Q of a 3-component spin field (j)(x) is the number of times (j)( x ) winds 
the sphere S 2 . It can be expressed as the integral over the space-time of a local operator Q(x): 

Q(x) = —e IMV e i jk4>i(x)du,<t>j(x)d v (f)k(x) ; (2) 
Q(x) is the divergence of a topological current |25|[p6|, 

Q(x) = d^K^x) . (3) 

All classical solutions with non-trivial topology - the fc-instantons - have been explicitly 
found p7||. At a quantum level, the only available prediction comes from the semiclassical ap- 
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proximation; the 1-loop resultjj for the size distribution of instantons in the physical vacuum is 
(l/V)dN/dp oc 1/p. This expression gives a logarithmic divergence in the ultraviolet limit p — > 0. 
Moreover, numerical approaches for the evaluation of the instanton size distribution fig| [^of indi- 
cate a growing behavior when the size decreases (the lowest size investigated is p ~ 0.2 £q). 

The topological susceptibility is defined as the correlation at zero momentum of two topological 
charge density operators Q(x) : 

X = J d 2 x (Q|T[Q(x)Q(0)]|0) . (4) 

The prescription defining the product of operators in Eq. (^) is |2£| 

(0| T[ Q(x)Q(0) ] |0) = 0„<O| T[ K^(x)Q(0) ] |0> . (5) 

This prescription eliminates the contribution of possible contact terms (i.e. terms proportional to 
the (5 function or its derivatives) when x — > 0. 



3 Scaling: mass gap measurements 

The general parametrization of the FP action of the 0(3) a-model is 

A FP {4>) = (5 { - \ J2 P(r)(l - m ■ cj>(x + r)) (6) 

x.r 

+ ^2 c(x 1 ,x 2 ,x 3 ,x 4: )(l - (f>(xi) ■ (j>(x 2 ))(l - (f>(x 3 ) ■ <j){Xi)) + ... }, 

where p represents the perfect discretization of the Laplacian and all the (infinite) multi-spin 
couplings with more than 4 fields have been only implicitly indicated. In this paper we test the 
properties of the 24-couplings parametrization reported in Table 4 of Ref . , which contains only 
one-plaquette terms with no more than 8 fields. 

This parametrized FP action should reproduce the continuum even at moderate values of the 
correlation length: cut-off effects on the spectral properties of the lattice theory should be absent, 
or at least strongly suppressed. This is confirmed by numerical evidences ||. We want to check 
this point further. 

We have performed our test by comparing alternative lattice definitions of mass gap. The 

3 Using the perturbative RG one can prove that this result is in fact asymptotically exact as p — > [H . 
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standard definition is related with the large-distance behavior of the wall-wall correlation function^ 

G w (y-x) = - ^2 G ( x ux;Vi>v) > ( 7 ) 



xi.yi 



where 



G(x) = {cj>{x) ■ 0(0)) . (8) 
The expected large-distance behavior, including periodic boundary condition effects, is: 



G w {x) 



A«, 



exp 



exp 



L — x 



for — > x > i u 



(9) 



It is also possible to define the diagonal wall-wall correlation function 



Get 



y-x 

v/2 



V2 



^2 G(x 1 ,x-x 1 ;y 1 ,y-y 1 ) 



whose large-distance behavior is 
A d 



G d (x) 



exp 



exp 



L/V2- 



for 



L 
2V2 



> x » £ d 



(10) 



(11) 



One of the scaling tests consists in verifying that the ratios £, d /£,w and A d /A w keep constant when (3 
varies. In addition, the comparison between £, d (A d ) and (, W (A W ) allows to test directly the rotation 
invariance: scaling violations are revealed by deviations from 1 of the ratios £, d /t; w and A d /A w . 

An alternative definition of the correlation length comes from considering the second moment 
of the correlation function 



e G 



We use the following lattice definition of £g : 

e 1 



Jd 2 x\x 2 G(x) 
Jd 2 xG(x) 

' 6(0,0) 

G(o, 1) 



4sin^(7r/L) 

where G(k) is the Fourier transform of G(x), given by 



(12) 



(13) 



(14) 



The zero component of G(k) is by definition the magnetic susceptibility \ m . In the scaling region 
the ratio (,g/(,w must be a constant, scale-independent number. The ratio should be equal to 1 
within 1% M. 



4 A small cut-off effect survives in the following definitions of the lattice mass gap since the ordinary lattice field 
instead of the FP field operator if> FP is used |2||5 1 
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Let's introduce the quantity Ac — Xm£,Q 2 £,w- The adimensional ratio A w /Ag is scheme- 
independent in the scaling region. We expect A w /Ag ~ 1, because the two-point functions should 
be almost saturated by the lowest energy state. 

We performed Monte Carlo (MC) simulations for several values of j3 corresponding to correlation 
lengths £ varying from ~ 2 to ~ 34 lattice units. We adopted a 4-hits Metropolis algorithm which 
is the only local algorithm available for actions with multi-spin interactions. The local nature of 
the algorithm will be essential in the application of the heating techniques to be discussed later. 
The acceptance of the pseudo-random changes of the Monte Carlo has been adjusted at 50%. 

In Table I we report values of correlation length and correlation function coefficient obtained, 
for the three definitions, at various /3's. The fits to G w and Gd have been performed choosing 
Xmin — we have checked that larger values of x m i n give results consistent within errors. Table II 
shows the ratios of these different definitions, analyzed using jackknife techniques. We can observe 
nice scaling and rotation invariance: in the region of /3 between 0.8 and 1.0 (£ between ~ 5 and 
~ 12) the uncertainty is lower than 1%. Moreover, with the same precision the ratios £,g/£,w and 
Aq/A w are consistent with one, so confirming the picture of an 0(3) model with a one-particle 
spectrum (see the discussion in Section 2). 

Table II shows also a test of asymptotic scaling, reporting the quantity (£g/(/?)2z) _1 , where 
f(P)2i — 27r/3e _2,r ^ is the 2-loop renormalization group function. In the asymptotic scaling region 
this quantity should be constant. However, numerical outcomes exhibit a steady growth and, 
moreover, they are far from the theoretical value Ma/App = 9.802494 (App is the A parameter 
of the theory regularized with the FP action) coming from the exact result for Mq of Ref. f3(i| 
and from the relation App = 8.17 A$t [||- We stress that this is not in contradiction with the 
properties of the FP action which do manifest only in the physical scaling of the lattice theory. Our 
hypothesis is that the deviation from asymptotic scaling is due to large non-universal corrections 
in the function f((3). We have fitted data in the seventh column of Table II taking into account the 
first non-universal correction: fsi(fi) — f(f3)2i(l + S/(2-Kf3j). Imposing for Ma/App the previously 
quoted theoretical value, we find 5 = —1.746(23), xVd.o.f. ~ 0.8. The higher order corrections 
to f (I3)2i are non-negligible in the region of /3 accessible to Monte Carlo simulations and it is 
essential to include them when the scaling of physical quantities is studied with this kind of action. 
Involving in the fit procedure also the ratio Mq/ App, we find Mq) 'App — 9.97(5), <5 = —1.80(3), 
X 2 /d.o.f. ~ 0.9. Nevertheless, we observe that the uncertainty of this measure is much larger than 
the naive error of the fit here quoted. A more reliable error can be found by checking the stability 
of the result under various alternative fits. Our conclusion is that the real uncertainty is about 
10%, being therefore our determination of Mq/ App compatible with the prevision of Ref. |$0| . 
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4 Stability of lattice classical configurations 



When instantons are discretized on the lattice, generally their action is no more scale invariant. In 
the case of the standard lattice action, for instance, the action of instantons decreases with their 
size [^j and, as a consequence, instantons are not stable under local minimization of the action, 
but shrink up to destruction ||. 

The FP action is scale invariant (up to a minimum size (l3)) and so its topological classical 
solutions are stable under minimization. When a parametrized form of the FP action is used 
in numerical simulations the statement about stability of instantons becomes likely approximate: 
however, a huge improvement is expected with respect to the standard action, and with respect to 
the tree level Symanzik improved action Q too, where only 0(a 2 ) cut-off effects are removed. 

Anyway, we did not study the absolute stability of instantons under minimization, our main 
interest being to check to what extent a fixed (finite) number of steps of local minimization modifies 
the lattice topological structures. This is a way to have control over the loss of topological signal 
which affects the lattice determination of the topological susceptibility by cooling techniques 
(see Section 7). 

In view of this, we have discretized on a 60 2 lattice the exact 2-instanton continuum solution 
defined on a torus p8|p| (1-instantons do not exist on a torus [Q(l^]). The study consisted in 
the observation of the behavior during 100 steps of cooling of both the topological charge and 
the action of the discretized 2-instanton. We studied 2-instantons of sizeP] ranging from 0.2a to 
5a, comparing the performances of the FP and tree level Symanzik improved actions. We have 
adopted an improved version of the lattice charge operator (see Section 7). The local minimization 
of the lattice action was obtained by making use of the Metropolis algorithm in the usual fashion, 
but with the request that only upgradings which lower the action are accepted. A sweep of this 
P = oo Metropolis is a cooling step. As a preliminary task we have checked that the efficiency of 
the algorithm of minimization is independent of the action; this test has been realized by cooling 
lattice configurations with trivial topology, and verifying that the behavior of charge and action is 
similar in the two cases. 

The results of the study on the 2-instanton are summarized in Figs. 1 and 2. We observe that in 
the case of the FP action topological charge loss is tiny if the size of the 2-instanton is larger than 
0.9a. For size smaller than 0.9a the discretized configuration becomes unstable and it is destroyed 
by few cooling steps. The deviation from 2 of the topological charge of the cooled instantons with 
size smaller than 3a is explained by a residual scale dependence of our (only partially improved) 

5 Wc thank M. Blatter and R. Burkhalter for many helpful suggestions in this regard. 
6 The size is defined as in the case of instantons on the sphere 



G 



definition of the lattice charge operator. In the case of the Symanzik action qualitatively the same 
behavior is observed, with the important difference that the topological signal is almost totally 
lost already at size 1.5a. 

A detailed analysis of the behavior of the charge during the cooling procedure gives some indi- 
cations about shrinking effects. Shrinking can be observed in both the Symanzik and parametrized 
FP actions, increasing at small sizes (where the scale symmetry breaking of the discretization is 
amplified); however, it is much reduced in the case of the FP action (see Table III): at p = 2a, for 
instance, it is ~ 2 orders of magnitude smaller. 

In Fig. 3 we compare the different cooling shapes of a small-size 2-instanton (p = 1.2a) when 
the FP, the standard and the Symanzik actions are used; being the algorithm of minimization non- 
deterministic, we have checked to what extent results fluctuate under changing of the sequences of 
random numbers used in the minimization algorithm. 

Although single-instanton classical solutions are not allowed on a torus, it is nevertheless inter- 
esting to investigate the amount of topological loss under cooling in the Q — 1 topological charge 
sector, which is dominant among classical configurations at thermal equilibrium. It is impossible 
to perform a study strictly equivalent to the previous one in the Q=2 sector, because 1-instanton 
solutions are missing on a torus: in their place, we have used 1-instanton solutions with constant 
boundary conditions (^?J which are quasi-solutions for p/L <C 1 (L is the lattice size). We have 
found substantially the same results of Figs. 1 and 2 as far as the minimum sizes are concerned: 
the main difference is that the observed action is slightly above the expected value Air at large p, 
going asymptotically to 4tt when p/L decreases (in Fig. 4 we report the results for the FP action); 
clearly this is a boundary effect (the configuration does not satisfy the correct periodic boundary 
conditions). 

5 Topology: the geometrical method 

In this Section we check how the geometrical method for the determination of the topological 
susceptibility Q works when the FP action is used. 

The geometrical definition Q g of the lattice topological charge is an attempt to recover the 
topological structure of the continuum theory by interpolating the discrete lattice field with a 
continuous field on which the topological charge can be measured. In formulae: 



where x* indicates a vertex of the dual lattice (corresponding to a plaquette of the 2-d lattice) 




(15) 
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and 

= 4^ {MX^i, 4»M + {<jA){4>ufo,<t>i)} . (16) 

4>i, 4>2, 03, 4>a are the four spins belonging to the plaquette x* and (aA)(<f>i, fa, 4>z) denotes the 
signed area of the spherical triangle cf>i , 4>2 , <fe Q 

The topological susceptibility constructed in terms of the geometrical charge Xg — (Q^/V) @ 
fails to reproduce the correct continuum behavior: the reason is that the geometrical definition as- 
signs non-zero topological charge to lattice structures of unphysical size (dislocations) . In the case 
of the standard lattice action, dislocations proliferate in the f3 — > oo limit |||)| since their action is 
definitely lower than 47r, the minimal (physical) value in the Q — 1 sector. As a consequence, they 
produce a non-scaling topological signal. Now, it is interesting to investigate to what extent the 
situation improves with the FP action, checking, in particular, if dislocations are still present and, 
if so, if their action is at least closer to 4ir. In order to have an indication about this subject, we 
have studied the cooling of a discretized 1-instanton solution (with constant boundary solutions 
to infinity). This lattice configuration - which is not a local minimum of the action on a toroidal 
lattice - is driven by the cooling algorithm to the nearest minimum in the Q g = 1 sector (lying on 
the border between the Q g = 1 to the Q g = sectors), and eventually it passes into the Q g = 
sector. The action of this border configuration is, in the case at hand (see Fig. 5), S/An = 0.92 < 1, 
thus indicating a dislocation. This outcome (though the method was not rigorous) suggests that 
the action of dislocations is closer to Air compared to the case of the standard action, where the 
minimal action in the Q g = 1 sector is 0.53 H. Anyway, dislocations are still present, so, apart 
from non-scaling effects due to the UV topological dominance of the model [|l9| |^(| , a deviation of 
Xg from the renormalization group behavior is expected. 

We have performed numerical simulations on Xg ( see Fig. 6 and Table IV) finding Xg ~ 
/3 2 e _47rQ/3 , with a = 0.8315(11), x 2 /d.o.f. ~ 0.8. For comparison, we report the a coefficient 
obtained in the case of the standard action |7|, a = 0.625: clearly the effect of dislocations is lower 
with the FP action. 



6 Topology: the field theoretical method 



7 One can demonstrate that 

exp Qiovt^ = p _1 {1 + <f>i ■ <f>2 + <j>2 ■ <i>3 + 03 • 4>i + ■ (<t>2 x 03)} , 

p= {2(1+01 -<t>2) (1 + 02 -<f>3) (1 + 03 ■0l)} 1/2 > , (17) 
thus providing an explicit formula for Q g . 
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6.1 Perturbative calculations 



We follow the field theoretical method |15| for the determination of the topological susceptibility 
from the lattice. First of all, a lattice topological charge density operator is defined as a local 
operator having the appropriate classical continuum limit |3l| ; our choice is |0 ] : 

Q L ( X ) = -^- e ^ijk4>i(x) (4>j( x + M) - <Pj( x - ^)){^k{x + v) - (j) k {x-v)) . (18) 
oZn 

From the prescriptions of field theory it comes that a finite multiplicative renormalization 
connects the matrix elements of the lattice topological charge density with those of its continuum 
counterpart: 

Q L {x) = a 2 Z(/3)Q(x) + 0(a 4 ) . (19) 
We recall: /3 = 1/g. The lattice-regularized version of \ is 

= ( £ Q L WQ L (o) ) = 4 ( ( £ q l w n > ( 2 °) 



x L 



where L is the lattice size. 



A prescription equivalent to Eq. (|j) does not exist on the lattice and therefore the contribution 
of the contact terms must be isolated and subtracted. These contact terms appear as mixings with 
the action density S(x) and with the unity operator I, which are the only available operators with 
equal dimension or lower. In formulae 

X L (P) = a 2 Z{pf X + a 2 A(fi) (S(x)) + P(J3) (I) + 0(a 4 ) . (21) 

In Eq. ( pl| ) the quantity (S(x)) is intended to be the non-perturbative part of the expectation 
value of the action density, i.e. it is a signal of dimension two. 

The extraction of the physical value of the topological susceptibility from numerical simulations 
requires the evaluation of the renormalization constants Z((3), A(j3), P{fi). 

The standard way to proceed is to perform the calculation in perturbation theory. We have 
decided to consider in our calculations the quartic terms of the parametrization of the FP action 
with coefficient 0(1CP 2 ) or larger. The procedure followed for the calculation of Z{j3) and P((3) is 
illustrated in Ref. ||; for A{0) we refer to Ref. @. 

An intermediate step in the calculation of Z{0) is the evaluation of the renormalization con- 
stants of the fields ir(x) and of the coupling constant g. We have found them to 1-loop order by 
imposing 

Z^(g,tm)T^(g,h,a;p) = T^(g r ,h r ,^p) , (22) 
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where g = Zg IS g r , h = Zg' IS Z^ IS h r . We quote here the results: 

Zf%,iua) = 1 - 2L((ia)g + 0(g 2 ) , (23) 
Zjs{g,na) = 1 - (L(fxa) + Cl )g + 0(g 2 ) , 

where L(x) — hxx + In 2 — y, being y = 0.0382 and c\ = 0.1262. We observe that at this 
order Z^f s is independent from the couplings of the multi-spin terms which affect only c\. 

For the calculation of Z(/3) we compare the 2-point T-functions with one operator insertion on 
the lattice and on the continuum: 

Z(fl- 1 z]™(0,»a)rZ imr {g,h,a;p,q) = T^(g r ,h r , mP ,q) . (24) 

^Q,-mr does not take contribution from the couplings of multi-spin terms at 1-loop order, so our 
1-loop determination of Z(f3) is exact. We show the result: 

Z(/3) = l + ^+0((3~ 2 ) , z x = -0.94237 . (25) 

P 

We calculated the perturbative tail P(/3) = J2 n =iPn/ P n U P to four loops. The results for 
infinite volume are: 

P 4 = 1.97429 x 10~ 4 , p 5 = -2.24879 x 10~ 4 . (26) 

Only is affected by our approximation. These numbers were extracted by performing numerical 
integrations at finite volume and extrapolating the results to infinite volume. 

We have calculated the 03 coefficient in the perturbative series A(j3) = J2 n =3 a n/ P n by com- 
paring r(Q£)2/y !ir7r with Ts^tt at the lowest order (for details see the Appendix in Ref. p0|]). The 
value of 03 extrapolated to infinite volume is 7.76206 x 10 -4 . This result also is independent from 
the multi-spin terms. As in the case of the tree level Symanzik improved action [pof , the mixing 
of the topological susceptibility with the action density is absolutely irrelevant. 

6.2 Numerical simulations 

We have performed Monte Carlo simulations for \ L on a 60 2 lattice over an extended region of 
values of f3, the algorithm being the same 4-hits Metropolis of previous simulations. In Table V 
we report the results: binning techniques have been applied in order to take account of the effect 
of correlations in the evaluation of errors. In order to control finite size effects, we have repeated 
the simulations at = 1.0, 1.1, 1.2 with the care that L/£g ~ 6. We observe (see Table V) that 
only the datum at = 1.2 appreciably differs from the result on the 60 2 lattice. 
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In the range of values of (3 accessible to our lattices, the perturbative determinations of the 
previous subsection are not sufficient to have a reliable estimate of the renormalizations Z{0) and 
P{f3)- Moreover, the standard technique to estimate some further terms in the perturbative expan- 
sion of P(f3) by a fit on numerical data at (3 higher than a certain f3 t (the so-called perturbative 
region) does not work: the extrapolation of P((3) to lower values of /3, where the physical signal 
is still living, is not reliable, since this procedure is unstable under changes of the conditions of 
the fit. So, we followed an alternative method for the determination of the renormalizations of 
the lattice topological susceptibility: the "heating" method J3^|33]]. The heating method, rely- 
ing only on MC techniques, is fully non-perturbative. It consists (for details see also |L2]] |2(]] ) in 
constructing on the lattice ensembles of configurations {Ct}, each configuration of the ensemble 
being obtained by performing a sequence of t local Monte Carlo sweeps starting from a discretized 
classical configuration Co - a large instanton or the flat configuration. For small t, the heating 
process thermalizes only the small-range fluctuations which are responsible for the renormaliza- 
tions: when the starting configuration is a large instanton (flat configuration), measuring Q L (x L ) 
on the ensembles {Ct}, a plateau at the value of Z({3) (P(f3)) is observed after a certain time, not 
depending on /?, corresponding to the thermalization of quantum fluctuations. 

In Table VI we show the values of Z(JS) evaluated by the heating procedure at various /3. 
Data reported therein have been obtained by heating a discretized 1-instanton with constant 
boundary conditions^] (size 10 lattice units) on a 60 2 lattice. The ensembles {Ct} have been 
obtained performing t sweeps of the Metropolis algorithm. In Table VI we also compare the 
non-perturbative determinations of Z{f3) with our 1-loop perturbative calculation: the latter well 
reproduces the non-perturbative results at the largest values of /3, but a deviation is observed at 
the smallest ones, where the 1-loop perturbative expansion is expected to fail. Attributing this 
discrepancy to the next -to-leading terms of the perturbative expansion of Z((3), we have found by 
a fit on data for j3 < 1.4: z 2 = 0.473(7), z 3 = -0.234(9) (x 2 /d.o.f. ~ 0.1). 

In Table VII we show the outcomes of the heating method for P(f3) at various (3. Data have 
been obtained on a 60 2 lattice. Using data from heating, we have fitted two further coefficients, 
Pq and p7, of the perturbative expansion of P(f3); the result is: pa — (7.54 ± 0.09) x 10~ 5 , pj = 
(3.02 ± 0.10) x 10 -5 , with x 2 /d.o.f. ~ 0.9. In Fig. 7 we compare three different determinations 
of P((3): the first is the pure 4-loop perturbative calculation; the second is a fit which combines 
the perturbative calculation with the equilibrium values of \ L & t values of (3 larger than 1.6; the 
last is the non-perturbative determination of P{/3). As it can be clearly observed, the curve from 
the second determination badly reproduces data from the heating method for the values of (3 not 
included in the fit, i.e. (3 < 1.6. This explains why the standard procedure of extracting the 
perturbative tail joining perturbative calculations and MC data from the perturbative region fails 
in this case. 

8 For the purposes of this method it is sufficient to thermalize around any smooth configuration with non-trivial 
topology. 
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Having obtained a reliable determination of the renormalizations Z(/3) and P(/3), we are now 
in a position to extract from \ L the physical quantity x ( see Eq. - the term proportional to 
A{(3) is negligible). 

In Table VIII (see also Fig. 8) we report, for different values of f3, the quantities x/^fp an d 
x/Mq (in the latter entry we use the relation between Kpp and Mq which can be deduced from 
Ref. we also report our determination of £g, as obtained from the fit of Section 3 (the one 

with the ratio Mq/ Kpp fixed at the theoretical value of Ref. pc|). 

In Table IX we compare our present results for x/Mq with other determinations known in 
literature. Data in columns (2) and (3) have been obtained in the field theoretical approach with 
the standard definition of the topological charge operator, but using other lattice regularizations 
of the action [[l2| |34| . In column (4) a parametrized FP action and the FP lattice charge operator 
have been used. The FP charge operator has been constructed by the application of one RG 
transformation on the geometrical charge operator |l3| . A substantial agreement is observed among 
determinations using the standard charge operator in the field theoretical approach, while data 
for x/Mq of column (4), where a FP charge operator is used, are systematically larger. We argue 
that the disagreement originates from the different performance in the short distance regime of 
the two lattice definitions of the charge operator. Indeed, the field theoretical definition is likely 
to underestimate the topological content of small size topological configurations (as Fig. 1 clearly 
singles out), while such loss of topological signal is definitely less relevant with the FP operator of 
Ref. Jll| (Fig. 2 of the previously quoted reference) . This different behavior strongly reflects upon 
the determination of the topological susceptibility, owing to the UV topological dominance of the 
model. 



7 Topology: the cooling method 

In this Section we turn to an alternative way to face the matter of renormalizations, known in 
literature as the cooling method || [|l6|] . The method relies on a local algorithm of minimization of 
the action of lattice configurations - the cooling; the purpose is to destroy the quantum fluctuations 
of the lattice configuration generated by the MC algorithm (which give rise to the renormalizations 
Z{(5) and P(/3) of Eq. (|2l|)), trying to preserve the background topological structure. If the 
background structure is an extended classical configuration, after a large enough number of cooling 
iterations, a measure of the lattice topological charge on the cooled configuration well approximates 
an integer^ 

9 If the background configuration is not slowly varying, scale effects produce deviations from integer values, due 
to the scale dependence of the field theoretical charge operator (see Section 4). 
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As we have seen in Section 4, shrinking effects are tiny with the (parametrized) FP action, so 
the most important drawback of the cooling procedure is overcome. 

In Table V we report the results of the simulations on a 60 2 lattice. The minimization algorithm 
has been described in Section 4; the "cooled" equilibrium value of the topological susceptibility 
has been determined by measuring the topological charge on configurations obtained from the 
equilibrium ones after 30 iterations of the algorithm. 

Here we adopt a different definition of the lattice charge operator. It consists in a tree level 
Symanzik improved operator, where the first three irrelevant terms have been eliminated: 

Q L ( X ) = -^-£v e ijk(t)i(x)D tJ ,(j)j(x)D l ;(j) k (x) ; (27) 
is the improved lattice derivative: 

n _ 1225 wt 1 ) _ 245 vy(3) , 49 tt( 5 ) 5 nP) /oo\ 

A* — 1024 V ^ 1024 V ^ ~ t ~ 1024 V A' 1024 V M ' \ A0 > 

Vji n) = ^(4>i(x + nfx) - (pi(x - nfj,)) . 

Being the first irrelevant terms absent, this operator is less sensitive to quantum fluctuations than 
the definition of Eq. ( |l8| ) and, as a consequence, the determination of the topological charge of the 
cooled configurations improves. 

For each value of /3 typically 100000 configurations were generated, the cooling procedure being 
performed once every 100 configurations. In Table VIII and in Fig. 8 we compare the cooling 
determinations with the results of the field theoretical method. Data from cooling nicely reproduce 
the outcomes of the field theoretical method, so confirming that topological structures are affected 
in a negligible way by the cooling procedure. 

The results for x/Mq obtained on large lattices by cooling (see Table VIII) allow now to 
appreciate a slight deviation from scaling - a non-scaling behavior is observed also in Ref. 
Such deviation is expected since a (/3-dependent) part of the whole topological signal is lost on 
the lattice due to the UV dominance of the model. However, the discrepancy of our data with the 
results of the previously quoted reference is still present and so it cannot be explained by finite 
size effects. 
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TABLE I: correlation length £ and correlation function coefficient A with the FP action. 



(3 


L 








A G 


A m 


A d 


0.6 


36 


2.14(4) 


2.188(15) 


2.180(21) 


5.15(18) 


4.84(4) 


4.84(6) 


0.7 


36 


3.14(5) 


3.16(5) 


3.17(5) 


6.3(1) 


6.22(11) 


6.19(12) 


0.75 


36 


3.89(7) 


3.89(7) 


3.89(6) 


7.04(7) 


6.96(12) 


7.00(9) 


0.8 


36 


4.81(7) 


4.83(9) 


4.80(6) 


8.103(23) 


8.01(14) 


8.11(8) 


0.85 


36 


5.87(9) 


5.86(13) 


5.97(9) 


9.37(3) 


9.45(17) 


9.30(10) 


0.9 


48 


7.34(16) 


7.34(21) 


7.34(15) 


10.98(3) 


10.99(25) 


11.12(14) 


1.0 


48 


11.71(29) 


11.8(4) 


11.89(27) 


15.01(14) 


15.0(3) 


15.4(3) 


1.1 


130 


18.5(2.1) 


18.1(2.7) 


19.1(2.1) 


22.3(4) 


24(4) 


22.2(1.7) 


1.15 


150 


24.4(1.9) 


24.3(2.5) 


23.7(1.6) 


26.5(3) 


26.7(2.5) 


28.2(1.2) 


1.2 


150 


34(3) 


34(4) 


33.0(2.9) 


33.1(1.3) 


32(3) 


35.1(2.6) 



TABLE II: ratios of different definitions of correlation length and correlation function coefficient; 
in the last column we show (£g/(/?)2/) _1 - 



a 


L 






A W /A G 


Ad/ A w 


(e G /(/3)2i)- 1 


0.6 


36 


0.976(18) 


0.997(7) 


0.94(3) 


1.001(11) 


5.59(23) 


0.7 


36 


0.994(5) 


1.003(9) 


0.987(11) 


0.996(17) 


5.88(10) 


0.75 


36 


0.9977(21) 


0.998(11) 


0.988(9) 


1.006(18) 


6.08(11) 


0.8 


36 


0.997(4) 


0.994(11) 


0.989(14) 


1.011(19) 


6.28(9) 


0.85 


36 


1.002(6) 


1.019(15) 


1.009(19) 


0.984(22) 


6.66(11) 


0.9 


48 


1.000(6) 


0.999(18) 


1.001(22) 


1.012(27) 


6.88(16) 


1.0 


48 


0.993(12) 


1.008(25) 


1.00(3) 


1.03(3) 


7.28(18) 


1.1 


130 


1.02(4) 


1.06(11) 


1.09(15) 


0.91(15) 


7.8(9) 


1.15 


150 


1.00(3) 


0.98(6) 


1.01(10) 


1.05(10) 


7.8(6) 


1.2 


150 


0.99(5) 


0.96(10) 


0.97(13) 


1.09(16) 


7.5(6) 
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TABLE III: topological charge (in the field theoretical definition) of 2-instantons of different initial 
size at the beginning and after 100 step of the cooling procedure with the Symanzik and the FP 
actions. The lattice size is 60. 





Qo 


QlOO 


p/a 




Symanzik action 


FP action 


4 


1.99856751 


1.99856600 


1.99856713 


2 


1.93896373 


1.92357210 


1.93858846 


1.2 


1.64901344 


0.00001479 


1.64991933 



TABLE IV: Xgi Xg/^FP an d Xg/^c f° r var i°us (3 and £g on a 60 2 lattice. £g is obtained from 
the fit of Section 3 (the one with the ratio Mq/Afp fixed at the theoretical value of Ref. IpOfl ). 



(3 


Cg 


Xg X 10 4 




XglM 2 G 


1.05 


15.42(7) 


8.27(26) 


18.9(8) 


0.197(8) 


1.1 


19.82(9) 


5.63(21) 


21.2(1.0) 


0.221(10) 


1.15 


25.58(11) 


3.61(9) 


22.7(8) 


0.236(8) 


1.2 


33.13(13) 


2.25(7) 


23.7(9) 


0.247(10) 
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TABLE V: \ L an d Xcooi versus (3 on a 60 2 lattice. Data with *, **, and *** arc obtained on lattices 
70 2 , 130 2 and 180 2 , respectively. 



(3 


X L x 10 5 


Xcool x 10 4 


(3 


X L x 10 5 


0.8 


27.0(4) 


- 


1.5 


1.91(6) 


0.85 


23.8(3) 


- 


1.52 


1.94(9) 


0.9 


20.30(29) 


- 


1.55 


1.77(8) 


0.95 


17.13(18) 


9.2(5) 


1.58 


1.61(7) 


1. 


14.12(16) 


6.42(29) 


1.6 


1.54(5) 


*1. 


13.97(20) 


5.98(27) 


1.65 


1.39(6) 


1.05 


11.28(12) 


4.48(21) 


1.7 


1.23(5) 


1.1 


9.18(10) 


2.79(19) 


1.75 


1.06(6) 


"1.1 


9.48(25) 


2.9(3) 


1.8 


0.98(5) 


1.12 


8.01(9) 


2.35(13) 


2. 


0.68(3) 


1.14 


7.36(9) 


1.88(11) 


2.2 


0.487(15) 


1.16 


6.72(8) 


1.48(9) 


2.3 


0.426(14) 


1.18 


6.05(8) 


1.33(10) 


2.5 


0.326(10) 


1.2 


5.42(6) 


1.00(9) 


2.8 


0.199(6)) 


^ 2 


6.75(25) 


1.27(13) 


3. 


0.160(5) 


1.22 


5.3(3) 




3.5 


0.095(3) 


1.225 


4.91(22) 




4. 


0.0576(9) 


1.25 


4.50(17) 




4.5 


0.0384(10) 


1.3 


3.82(15) 




5. 


0.0256(7) 


1.32 


3.35(10) 




6. 


0.0131(9) 


1.35 


3.06(15) 




7. 


0.00683(18) 


1.38 


2.96(15) 




8. 


0.00409(11) 


1.4 


2.72(23) 




9. 


0.00250(7) 


1.42 


2.40(7) 




10. 


0.00178(5) 


1.45 


2.36(18) 
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TABLE VI: Z((3) versus (3. Z(f3)u is the multiplicative renormalization calculated to one loop. 
Z(/3)mc is the multiplicative renormalization calculated by heating a discretized Q = 1 smooth 
configuration; the size of the lattice is 90 for (3 > 1.5 and 60 for the remainder; the statistic of the 
simulations is 1000. Since data on the plateau are correlated, we have estimated errors by jackknifc 
techniques. 



p 


Z{(3)u 


Z(P)mc 


5.0 


0.812 


0.8247(7) 


4.5 


0.791 


0.8096(8) 


4.0 


0.746 


0.7865(11) 


3.5 


0.731 


0.7599(13) 


3.0 


0.686 


0.7254(18) 


2.5 


0.623 


0.6804(25) 


2.1 


0.551 


0.631(3) 


1.9 


0.504 


0.591(4) 


1.7 


0.446 


0.566(4) 


1.5 


0.372 


0.514(6) 


1.4 


0.327 


0.479(9) 


1.35 


0.302 


0.466(10) 


1.3 


0.275 


0.450(10) 


1.25 


0.246 


0.430(11) 


1.2 


0.215 


0.406(12) 


1.15 


0.181 


0.380(13) 


1.1 


0.143 


0.362(14) 
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TABLE VII: P((3) np versus [3. P(/3) np is the mixing with the unity operator I calculated by heat- 
ing the flat configuration; the size of the lattice is 60, the statistic of the simulation is 1000. 



(3 


P(0) np x 10 5 


(3 


P{f3) np x 10 5 


1. 


7.58(24) 


1.25 


3.30(8) 


1.05 


6.48(19) 


1.3 


2.84(6) 


1.1 


5.49(14) 


1.35 


2.53(6) 


1.15 


4.60(11) 


1.4 


2.23(5) 


1.2 


3.84(9) 







TABLE VIII: x/^fp an d x/Mq for various (3 and £g for both the field theoretical and the cooling 
determinations. L is the size of the lattice. 



(3 


L 


£g 


Xl^FP 


XIMI 


Xcool/ A^p 


Xcool /M% 


0.8 


60 


4.74(3) 


8(7) 


0.09(7) 






0.85 


60 


5.92(4) 


9(4) 


0.09(4) 






0.9 


60 


7.46(4) 


8.9(2.5) 


0.093(27) 






0.95 


60 


9.45(5) 


9.4(1.9) 


0.098(20) 


7.9(5) 


0.083(5) 


1. 


60 


12.04(6) 


10.0(1.7) 


0.104(18) 


8.9(5) 


0.093(5) 


1. 


70 








8.3(5) 


0.087(5) 


1.05 


60 


15.42(7) 


10.3(1.5) 


0.107(16) 


10.2(6) 


0.106(6) 


1.1 


60 


19.82(9) 


11.3(1.5) 


0.118(15) 


10.5(8) 


0.110(8) 


1.1 


130 








11.1(1.3) 


0.115(13) 


1.12 


60 


21.94(10) 


10.3(1.4) 


0.108(14) 


10.9(7) 


0.113(7) 


1.14 


60 


24.31(10) 


10.7(1.4) 


0.111(14) 


10.7(7) 


0.111(7) 


1.16 


60 


26.93(11) 


10.9(1.4) 


0.113(14) 


10.3(7) 


0.107(7) 


1.18 


60 


29.86(12) 


10.6(1.4) 


0.111(15) 


11.4(1.0) 


0.118(10) 


1.2 


60 


33.13(13) 


10.0(1.3) 


0.104(14) 


10.5(1.0) 


0.110(10) 


1.2 


180 








13.4(1.5) 


0.139(15) 
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TABLE IX: Comparison of different determinations of x " £g = x/Mq- theoretical method 

and FP action (column (1), present work), field theoretical method and Symanzik action (column 
(2), fl2| ), field theoretical method and CP 1 standard action with explicit gauge degrees of freedom 
(column (3), p4|), FP charge operator and parametrized FP action^ (column (4), Jl3[). 



£g 


(1) 


(2) 


(3) 


(4) 


4.74(3) 


0.09(7) 


- 




- 


5.92(4) 


0.09(4) 






- 


6.057(17) 


- 


- 


- 


0.1004(9) 


7.46(4) 


0.093(27) 








9.45(5) 


0.098(20) 


- 


- 


- 


12.04(6) 


0.104(18) 


- 




- 


12.16(3) 




- 




0.1448(15) 


15.42(7) 


0.107(16) 


- 




- 


15.9(9) 






0.11(3) 




18.56(8) 




0.117(5) 






19.82(9) 


0.118(15) 








20.40(9) 








0.1893(27) 


21.2(1.1) 






0.12(3) 




21.94(10) 


0.108(14) 








21.96(9) 




0.121(7) 






24.31(10) 


0.111(14) 








24.57(10) 




0.124(7) 






26.93(11) 


0.113(14) 








29.10(13) 




0.126(7) 






29.86(12) 


0.111(15) 








32.58(13) 




0.127(8) 






33.13(13) 


0.104(14) 








34.4(3) 








0.224(5) 


38.60(16) 




0.121(10) 







We are indebted to M. Blatter and R. Burkhalter for the data relative to Fig. 2 in ll. 
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FIGURES 



FIG. 1. Cooling of small-size 2-instantons with the FP action. Action (dotted line) and 
topological charge (solid line) of cooled discretized 2-instantons are represented versus their size 
in lattice units. The cooled configurations have been obtained by performing 100 cooling steps on 
the discretized instanton solution. The action is reported in units of 4ir. The lattice size is 60. 

FIG. 2. As in Fig. 1 with the tree level Symanzik improved action. 

FIG. 3. Stability under cooling of a small-size 2-instanton. The topological charge of a cooled 
discretized 2-instanton is represented versus the step of cooling when the FP action (circle), the 
tree level Symanzik improved action (diamond) and the standard action (triangle) is minimized. 
The size of the 2-instanton in lattice units is 1.2. The lattice size is 60. The procedure has been 
repeated 100 times with different random sequences. 

FIG. 4. As in Fig. 1 with 1-instantons (constant boundary conditions). 

FIG. 5. Transition from Q g = 1 to Q g = during the cooling of a 1 instanton (constant 
boundary conditions) with the FP action. The behavior of action and topological charge in the 
geometrical definition is represented versus the cooling step. The instanton size in lattice units is 
15, the lattice size is 60. 

FIG. 6. Scaling behavior of the lattice topological susceptibility in the geometrical approach 
X g with the FP action. Diamonds represent data from MC simulations; the dotted line is the best 
fit, corresponding to an exponential decay with a ~ 0.83 (see the text); the solid line is the scaling 
behavior predicted by the renormalization group. The lattice size is 60. 

FIG. 7. Comparison between the perturbative and non-perturbative determination of P - the 
mixing of \ L with the unity operator - versus (3. The diamonds represent the MC data of x L , the 
circles the outcomes of the heating technique (the solid line is the relative fit); the dotted line is 
the 4-loop perturbative determination of P, while the dashed line is the best fit of the MC data 
in the perturbative region (/? > 1.6). 

FIG. 8. x/Mq versus £g, the correlation length in lattice units, with the FP action. 

(a) {x L {P) - P(P))Z((3)- 2 Z G (P) 2 on a 60 2 lattice (diamond); 

(b) x/Mq extracted by cooling on a 60 2 lattice (square); 

(c) the same quantity of (b) computed on lattices of different sizes keeping L /£g ~ 6 (filled square) . 
For the sake of legibility, data (b) have been slightly shifted in £g • 
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